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ON LINEAR DEGENERATE ELLIPTIC PDE SYSTEMS WITH 
CONSTANT COEFFICIENTS 

NIKOS KATZOURAKIS 


Abstract. Let A be a symmetric convex quadratic form on and U C R" 
a bounded convex domain. We consider the problem of existence of solutions 
u : fl G R" —> R^ to the problem 
N n 

(1) = fa, in fl, u = 0, on an, 

/3=1 = l 

when / S 1/^(0,R^). (1) is degenerate elliptic and it has not been considered 
before without the assumption of strict rank-one convexity. In general, it may 
not have even distributional solutions. By introducing an extension of distri¬ 
butions adapted to (1), we prove existence, partial regularity and by imposing 
an extra condition uniqueness as well. The satisfaction of the boundary condi¬ 
tion is also an issue due to the low regularity of the solution. The motivation 
to study (1) and the method of the proof arose from recent work of the author 
[K4] on generalised solutions for fully nonlinear systems. 


1. Introduction and the main result 

Let n, A > 1 be integers and consider a convex symmetric quadratic form A on 
the matrix space that is 

(1.1) ~ ^ ( ^ai[3j Qai QfHj — 0; Q € M 

Let also d M" be a bounded convex domain. We will follow the convention 
that Greek indices a,/3,7,... run in {1,...,A} and Latin indices run in 

{1,..., n}, even if their domain is not explicitly mentioned. In this note we consider 
the question of existence of solutions u : ft C M” —to the Dirichlet problem 

— fa: IH ft. Cl = 1, ..., A, 

/3dU 

u = 0, on dfl, 

when / G The operators DiUa, DfjUa will denote the i-th and ij-th 

partial derivatives of first and second order respectively of the Ua component of the 
map u = (ui,..., un)^■ 

It is well known that at least formally, (1.2) is the Euler-Lagrange equation of 
the convex functional 

(1.3) E{v,, n) .— / I o ^ ^ -^aifUj DiUa EjV,^ A ^ ( f'y U-y 


Key words and phrases. Degenerate elliptic 2nd order systems, Euler-Lagrange equation, Cal¬ 
culus of Variations, rank-one convexity, theory of Distributions. 
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in the Sobolev space i7g(r2,K'^). If A is strictly rank-one convex on that is 

when the following Legendre-Hadamard condition holds 

(1.4) ^ Aaifijriaaii^paj > c\r]\‘^\a\‘^, ry G a S K”, 


for some c > 0 (where | • | denotes the Euclidean norm on both and M”), then 
it is textbook material that the problem (1.2) has a unique weak solution which 
is a minimiser of (1.3) in the space Hq{Q,'R^), see e.g. [D, GM]. Moreover, by 
standard regularity results it follows that the solution is actually strong, lies in the 
space n iJg)(f2,K^) and satisfies the system a.e. on (e.g., [GM]). 

The primary advance in this paper is that we prove existence of solution to (1.2) 
without assuming the standard strict ellipticity condition (1.4), but instead only 
the degenerate ellipticity condition (1.1) and an extra constraint on the (nontriv¬ 
ial) nullspace of A which we explain later. Without (1.4) the functional (1.3) is 
convex but non-coercive and standard variational/PDE methods fail. This is not 
a technical weakness, since as we show by examples (Ex. 2, 3) the solution in this 
case does not exist as an element of any standard Sobolev space and may not exist 
not even in the sense of distributions unless a compatibility condition is satisfied. 
To the best of our knowledge, the problem (1.2) has not been considered before 
without the assumption of strict rank-one convexity. 

The idea of the proof is based on the vanishing viscosity approximation of (1.2) 
by the strictly elliptic systems 


(1.5) 


^ ( {Actipj -f — /q, 

P,i,j 


= 0 , 


in 11, a = 1,..., N, e > 0, 
on 911, 


and on the derivation of partial estimates along rank-one directions which are stable 
as e —^ 0. By introducing an appropriate variant of the Distributional solutions 
adapted to the degeneracy of A, we prove the existence of solution in this sense. 
We note that the satisfaction of the boundary condition is also a serious issue 
under the low regularity of the solution since the solution fails in general to be in 
Wjqj)( 11, and there is no general trace operator for ■^ioc(^>'^^) mappings. 

Before stating our existence result we need some preparation. Let A be a given 
tensor which satisfies (1.1) and will be fixed for the rest of the paper. The notation 

N{A : ^ 

will be used to denote the nullspace of A when A acts as a mapping 
3 Q ^ AQ := Y, {Ao.^pJQfiJ) e“ 0 e* G 


Evidently, {e*}, {e“} and {e“(8)e*} denote the Euclidean bases of K", R^ and R^^ 
respectively. Let us define the vector spaces 

n := A(A : ^ C 

E := span[{77 G : ?7 (g) a G II}] C R^. 


( 1 . 6 ) 
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The space 11 is the orthogonal complement of the nullspace of A (namely, the range) 
and contains the “rank-one directions of strict ellipticity”, that is 


^aiPoVaCiiilisaj > c| 77 p|ap, ? 7 (g)oenc 


D A/'n 


We will follow the convention that the same letters 11, S will denote the subspaces 
as well as the orthogonal projections on them. The meaning will be clear from 
the context, for example the projection map satisfies E = = E^ > 0 etc. The 

functional spaces of mappings / : C K" —E C valued in E will be denoted 
by E), W^’^(r2,E), C'°°(n,E), etc. Note also that we have 11'*- = {0} if and 

only if A defines a strictly convex quadratic form, whence we also have E = 
and n = in this case. Let now 

^(ff,E) := C'~(f],E) 

be the space of “test maps” valued in the subspace E C . We consider the space 
of Distributions “valued in E”, namely the dual space of E) 

^'(D,E) := (^(f^,E))*. 

We consider both spaces as being equipped with their usual topologies (which 

we will not need, so we refer to [F] for their definition). 

Definition. We will say that the map u : fl C M" —>• 
solution in ^'(fl, E) of the PDE system 


is a Distributional 




— faj IR 


when u S R-*^) and for all (j) € ^(fi, E) we have 


(1.7) 


/ 'y ^ Up — / y ( fa'P' 

r 2 _ ,o • • 




The following is our main result. 

Theorem 1 (Existence-Uniqueness-Partial regularity). Let n, A > 1 with Lt d K” 
a strictly convex bounded domain. Suppose that A is a quadratic form on 
which satisfies (1-1). We assume that the vector space 11 C R^" (i.e. the orthogonal 
complement of the nullspace of A, given by (1.6)^ is spanned by rank-one directions. 
Then, for any f G L^(fl, E), the problem 

y ^ -^aipj Ltj^jUp — fa, in D, O. — 1 ,..., A, 


tt = 0, 


has a unique Distributional solution u : LI C 


on. 


pN 


in ^'(D, E) which lies in 
L^(n, E). In addition, u is TL^~^\-dU-measurable on the boundary of any strictly 


Moreover, certain projections along rank-one directions of the Distributional gra¬ 
dient Du exist as functions: for any 77 0 a € H, the directional weak derivative 
of the projection Da{rj ■ u) exists in if (LI). 

In the above statement, denotes the (n — l)-Hausdorff measure, 

denotes the Euclidean inner product and “Dq” is the standard directional derivate 
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along a. We remark that / must be valued in the subspace S since this is a necessary 
compatibity condition arising from the degenerate nature of the PDE. 

The following considerations show that the results and the assumptions of The¬ 
orem 1 are optimal. 


Example 2 (Compatibity condition). The compatibility condition in Theorem 1 
which requires that f must be valued in E is necessary: the degenerate 2x2 system 

J Aui(a;i,a;2) = /i(a;i,X2), 

\ 0 = f2ixi,X2), 

has no solution whatsoever in any weak sense unless f 2 = 0. 


Example 3 (Partial regularity). In general the solution we obtain in Theorem 1 
can not be a Sobolev function. Let LI = Bi(0) C be the unit disc centred at 
the origin and choose a function f G which is not weakly differentiable with 

respect to xi for any X 2 . Then, the Dirichlet problem for the following degenerate 
elliptic single equation 

j Dl 2 u{xi,X 2 ) = f{xi,X 2 ), on LI, 

I It = 0, on dLl, 


has the solution 


/ f{Xi,S2)dS2dt2. 

-OO ^ —OO 


In the above, h is the function 


,, , (9{xi, ^Jl-xl)- g{xi,-sj\ - xj) \ 

:= -jx, 

g{xi, a/1 - xf) +g{xi,-^l - xf) 


where g G C^{dLl) is the function which is given by 

/ X2 pt2 

/ f{xi,S2)dS2dt2, xl+xl = l. 

-OO J —OO 

In view of our choice of f, the solution u is not in Wif^{Ll). 


Remark 4 (Nonuniqueness on the subspace of “degeneracies”). If we do not require 
the generalised solution we obtain in Theorem 1 to satisfy E"’~u = 0 then it may 
not be unique (unless A is strictly elliptic, in which case we have E-*- = {0}). Any 
extension of u from lf{Ll, E) to L'^{L1,M.^) which satisfies the boundary condition 
is also a solution. For instance, 

u := Eu -b E-^g, g G C°(Ll), 
is also a solution of the same Dirichlet problem for any g. 


The motivation to study the problem (1.2) and the method of proof come from 
the very recent paper of the author [K4] and its companion paper [K5]. In [K4] 
we proposed a new duality-free theory of generalised solutions which applies to 
fully nonlinear PDE systems. This approach allows for nonlinearities of any order 
and with discontinuous coefficients whilst the only a priori regularity requirement 
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of the solution is measurability. The standing idea of integration-by-parts which 
applies only to divergence systems is replaced by a probabilistic representation 
of derivatives which do not exist in the classical sense. Technically, this is done 
by utilising Young (parameterised) measures to describe the limiting behaviour of 
difference quotients over the compactification of the “state space”, that is the space 
wherein the derivatives are valued. 

Among other existence results in these papers, in [K4] we proved existence of a 
so-called “diffuse solution” to the problem 

^ g, f F{-,D^u) = /, in n, 

I u = 0, on dil, 


where F : Qx is a Caratheodory mapping and is the Euclidean 

space wherein the hessians D^u of smooth maps are valued. The result for (1.8) 
extends previous work of the author for (1.8) but on the realm of strictly elliptic 
systems and of strong solutions (see [Kl, K2, K3]). The crucial assumption for 
existence is a degenerate ellipticity condition which roughly requires F to be “close” 
to a linear degenerate system of the form we consider herein. This system has been 
solved in the “diffuse” sense in [K4] under assumptions stronger than those we 
consider herein, as a stepping stone in order to solve (1.8). To aim of this paper 
is to show that under the present weaker assumptions, the PDE system (1.2) has 
solutions in a certain distributional sense as well. 

A particular difficulty is the satisfaction of the boundary condition. In fact, 
the only reason that strict convexity of fl is needed is for the satisfaction of the 
boundary condition. The strictness is meant in the sense that dfl contains no non¬ 
trivial straight line segment. Although in general there is no trace operator because 
the solution may not be in any Sobolev space, yet it has “differentiable rank-one 
projections”. This means that the gradient Du does not exist as a whole, but only 
certain projections of it exist along rank one lines of Surprisingly, this suffices 

for a partial trace operator to exist. 

The following condition which was introduced in [K4] provides a sufficient con¬ 
dition about when the assumption of Theorem 1 that 11 is spanned by rank-one 
directions is satisfied (see (1.6)). 


Structural Hypothesis (SH) The tensor A satisfies (SH) if it can be decomposed 
as 

A . . _ d 1 /ll I I A-^ 

“r ’ ’ ’ “r 

and also: 

a) The symmetric matrices {A ^,..., A^} C M” are non-negative and the eigenspaces 

A^I^AT'- A^/: K" ^ K”) 

interseet for all 7 = 1,...,A^ along a eommon line. Here A)]) denotes the smallest 
positive eigenvalue of A'^. 

b) The symmetric matrices {B ^,..., B^} C are non-negative and have mutually 
orthogonal ranges. 


Note that (SH) trivialises when either = 1 or n = 1 since any symmetric non¬ 
negative matrix satisfies it. Although (SH) is quite restrictive, by its constructive 
nature is is evident how to demonstrate nontrivial examples when n, N > 2. 
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2. Proof of the main result 


Proof of Theorem 1. The proof is based on the approximation by the strictly 
elliptic systems (1.5) as e —)■ 0. 


Step 1. To begin with, fix / G S), e > 0 and consider the functional 


(2.1) E,iu,n) := 


E 


-^a.ij3j DjU^ 
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placed in Hq{Q,R^). Standard semicontinuity and regularity results (see e.g. 
[D, GM]) imply that the problem (1.5) has a unique strong solution u® G (i?^ n 
i?g(f2,K^) which also is a global minimiser of (2.1) over i/Q(r2,K'^). 


Step 2. Now we obtain estimates stable in e. 


We begin with some algebraic consequences of our assumptions. By (1.6), (1.1) 
and the Spectral theorem applied to the symmetric linear map A : —y K^", 

we have that there exists u > 0 such that 


(2.2) BAB = A = AB = BA, 

(2.3) ^ > u|BQ|2, QgK"^". 

a,i,P,3 

In (2.3), I • I denotes the Euclidean norm on namely IQp = J2a iQaiQai- 

The identities (2.2) say that A commutes with the projection map B on its range. 

(2.3) is a consequence of the fact that the restriction A|jj on the range is invertible. 
Moreover, we have that 

(2.4) XgK"^"' ^ G E C 

(2.4) says that the image of the linear map A : is a subspace of E. 

This is a result of the degeneracy of A and can be seen as folllows: every X G 
can be written as 

Nr? 

X = ^ (g) (g) G G K". 

A =1 


Then, by (2.2) we have (if the components of the projection B are denoted by Haipj 
and we use the symmetry of B) 


Nn^ 

Va ■— ^ail3j E! ^ai/3j ^/3 dj b 

0,ij A=l/3,i,j 

{ Nr? 


= E i E E df. 

K,k,i I A=1 P,j 


By the definition of E in (1.6), the above says that the vector rj := J2a dae“ belongs 
to E and hence (2.4) follows. 
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Now, since is a minimiser, by (2.2) and by the compatibility condition which 
says that S/ = / (because / is valued in E C R^), we have 


0 = > E,{u^,n) 


> 



ai^j E-iU^ EjUj^ 


j ^(S/)^u-. 


By (2.3) and by the symmetry of the projection map E, for any <5 > 0 small we 
deduce 


0 > 

> 



aijUj EiU^ EjUj^ 


f 5](e/),< 



E/7(S«^)7 


which gives 

(2.5) 0 ^ ^ 

^ Jn ° JQ. Jn 


Now we need a generalisation of the Poincare inequality and of the trace operator. 
The rest of the proof follows similar line to those of [K4], by we provide all the 
details for the convenience of the reader. 


Claim 5 (Partial Poincare inequality, cf. [K4]). Let n,n, E be as in Theorem 1. 
Then, there exists C = C{Ll,n, N) > 0 depending only on the diameter of fl and 
on the dimensions n,N such that, for any u G we have the estimate 

l|Sw||L2(n) < C'||nT>'“||L2(n)- 

We note that Claim 5 is actually true for any bounded open domain LI. 

Proof of Claim 5. Fix vectors e G and rj G and let us denote by the 

hyperplane normal to e. For any y G we set (see Figure 1) 

iy,e ._ g I g Qe |y g g-L I g ^ g ]^ • y _|_ ^g g 

We fix a function u G CQ(fl,R-^) and some x = y + te G L2. Then, for the projection 
T] ■ u along 77 we have 

\{'q ■ u)(y + te)\^ < f |e • D (?7 • u)( 2 /+ Ae)|^dA. 

By integration with respect to t € I*'’® and y gLE, Fubini’s theorem implies 

f \{r] ■ u){x)\'^dx < f f \e ■ D{r] ■ u){y + Xe)\'^d\\ d'H^~^{y) 

Jn Jq” \ Jiv.o ) 

< diam(n)^ / \rj ® e ■. Du{x)\^dx. 

Jn 


NIKOS KATZOURAKIS 



Since by assumption 11 C is spanned by rank-one directions of the form 77 0 e, 
by the definition of S in ( 1 . 6 ) the desired estimate follows by considering a basis 
of n consisting of rank-one matrices and by using a standard density argument in 
the Sobolev norm. □ 


Finally, by using Claim 5 and (2.5), by choosing i5 > 0 small enough we have the 
uniform in e estimate 


( 2 . 6 ) 


iSu®||L2(Q) -f < (^11/1^2(0) 


for all e > 0. By (2.6) and by weak compactness, there exist maps u S E) 

and U € L^(n,n) such that 


(2.7) 


Eu®—^ M, inL^(n,E), 
nuu®—^ c, mL^(n,n), 

along a sequence e = £fe—^Oasfc—>-00. 


Step 3. Now we prove existence of a solution to (1.2). Since for any e > 0 is 
a strong solution of (1.5), we have 

^ ^ ^ij^p ~ f cli 


a.e. on for all indices ol. Moreover, = 0 ■H"“^-a.e. on i9f2 by standard results 
on the pointwise properties of Sobolev functions and the regularity of the domain 
(see e.g. [EG]). Let us fix a test map S ^(fl, E). Integration by parts gives 

(2-8) [ ^crpj u^p D'fjK = f ~ ^ f 

By (2.4) we have that 

(2.9) ^ [ ^aiPj PgUp = 'y [ i^-^aiPj 

p,ij A,/3,z,j 

and hence by using that T,(j) = cj) and (2.9), (2.8) gives 

/ y ' {-^aiPj Elytra) = / y ( fjfpa ~ ^ I ^ ' (Eak^k) ^^X- 

Jn_ oJn Jn , 
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By letting s = Sk ^ 0, the convergences of (2.7) imply that the previously obtained 
limit map u G S) is a Distributional solution in ^'(D, E) of the system 

'y = fa. 

Step 4. We now consider the problem of the satisfaction of the boundary 
condition. The next result shows that for strictly convex domains we have a partial 
trace operator under the assumptions of Theorem 1. 

Claim 6 (Trace operator, cf. [K4]). Let D, If, E he as in Theorem 1. Then, there 
exists a closed -nullset E C dO, such that, for any T (s dO, \ E, we can find 
C = C{n,T) >0: 

l|Sw||L2(r.W"-i) < C* (l|Su||L2(n) + , 

for all maps u G In addition, we have that 

E = G dLl : [x + span[e\) H dQ = {a;}|. 

We note that Claim 6 is a minor extension of standard results (see e.e. [E, EG]). 

Proof of Claim 6 . Suppose E C dil is the closed set defined in the statement of 
the claim. By the strict convexity of D, it can be seen that TIT~^{E) = 0. Let us fix 
a function u G M.^) and a unit rank-one matrix 77 ® e S If C We cover 

9D \ E by a sequence of open cubes {Qj}^ whose sides are orientated parallel to 
{e, e-*-} (see Figure 2). For every cube Qj, we consider the sets 

Tij := Qj n n, Tj := Qj n on. 

Let us fix a triplet {Qj, Llj, Tj) and assume that e points towards the interior of Tlj 
(for otherwise we may replace it by —e). We may also restrict ry • u, D{rj ■ u) on Clj 
and define 

iv-u)j := {■n-u)xn,, {D{r] ■ u))j := D{iq-u)xn^. 

Then, standard estimates imply that for any x G Tj we have 

\{rj ■ u){x)\^ < C {j \{ri ■ u)j{xte)'^dt + J \(^D{r] ■ u))^{x + te) ■ e\'^dtj . 



Suppose now that Fj G W^’°°(]R") is a Lipschitz function such that when restricted 
to Dj, it satisfies the following: for each t > 0, the level set {Fj = t} coincides with 
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the translate of a portion of the boundary + te. Such a function Fj can be given 

by 

Fj{z) := sup{t > 0 : z € (fi + te) n z € 

and can be trivially extended to a Lipschitz function on K". Then, we integrate 
over X G Fj and apply the co-area formula (see e.g. [EG]) and Fubini’s theorem: 

\{V u)\^ ^ + \{D{ri-u)).{y) ■e\^^dn^-^iy)dt 

= C' / / (\ivu)j{y)\^+ \{D{rfu))dy)-e\^)dH'^-^iy)dt 

= C\DFj{x)\(^\{yu)j{x)\^ + \{D{yu)).{x)-e\^^dx 

< C'||I?Fj||i=o(R„) y (^|(?7 • ■u)(x)|^-b |? 7 G e : £)u(a;)|^) dx. 

By using the assumption that If has a basis consisting of rank-one matrices 77 G e 
and E is spanned by the respective directions rj G , the rest of the proof is an 
obvious application of a standard argument of partitions of unity. □ 

An application of Claim 6 shows that the Distributional solution u is 
measurable and u = 0 'H"“^-a.e. on dfl. In addition, it is ^-measurable on the 
boundary of any strictly convex subdomain of 

Step 5. We now show the uniqueness of the solution. Suppose the problem has 
2 solutions u, v and set w := u — v. Then, w G L^{Fl, E) is a solution in E) of 

N n 

Aaii 3 j = 0 , in D, w = 0 , on dVL. 

/3=ii.i=i 

(The satisfaction of the boundary condition is considered in the sense of the state¬ 
ment of the theorem.) Then, w can be extended on K" \ fl by zero as a function in 
L^(K”, E). By mollifying in the standard way, for any £ > 0 the mollified solution 
= w satisfies 

N n 

y] y] = 0 , onK", 

p=l i,j=\ 

and in® G E). Integration by parts and application of the inequality (2.3) 

give 

II IIl2(B") — ^ J Di{w )aDj{w )p = 0. 

Finally, by Claim 5 and by the compactness of the support of in® we obtain Ein® = 0 
on K". Since E-*-u; = 0, by letting e —)■ 0 we get in = 0 and as such we infer that 
the solution of the problem is unique. 

Step 6. In order to conclude it remains to show the weak differentiability of the 
projection 17 • u along the a-direction of K” when 17 C a £ 11. This is a consequence 


LINEAR DEGENERATE ELLIPTIC SYSTEMS WITH CONSTANT COEFFICIENTS 


11 


of the convergences in (2.7): for any (j) S we have 

{r]-u)Da(l)= lim {r] ■ u^'‘ )Da(t) 
Jn 


and hence 


= — lim 

k—>oo 


(l)Da{T] ■ 

= — lim / 6 (ti ® a : 

fc^oo ^ ’ 

/ (77 • u)Da,<^ = - / (f) [rj ® a ■. U), 

Jo, JQ 


for any (j) € As a consequence, we have that Da{r] -u) = r]^a : U and also 

that 77 0 a : [/ is in The proof of the theorem is completed. □ 
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